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(J> , We present our first attempts to formulate a thermodynamics-like description of 

CD 1 explosions. The motivation is partly a fundamental interest in non-equilibrium 

CD ' statistical physics, partly the resemblance of an explosion to the late stages of a 

^y\ ' heavy-ion collision. We perform numerical simulations on a microscopic model of 

^^ ' interacting billiard-ball like particles, and we analyse the results of such simulations 

f^^ , trying to identify collective variables describing the degree of equilibrium during 

f*^ , the explosion. 



1 Introduction 

73 . 

^ . The assumption of thermodynamic equilibrium at an intermediate stage of a 

heavy- ion collision is often incorporated in models of the colliding nuclear mat- 
ter. These models range from statistical models of nuclear multifragmcntation 
(>j " to the fluid dynamical models of the quark gluon plasma. In contrast, mi- 

'V^ . croscopic models of molecular dynamics type (e.g. RQMD, FMD and NMD), 

C^ ' which are based upon constituent interactions, do not contain this assump- 

tion. Such models are appropriate for testing to what extent thermodynamic 
equilibrium is actually achieved. And if it is not, the application of thermo- 
static concepts such as temperature and entropy becomes questionable. In 
this study we employ a very simple model, and focus on the thermodynamic 
or "overall" description of the system. 
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2 The model 

Our model consists of a number A of identical balls of radius r^c having mass 
m. They perform classical non-relativistic hard-sphere scatterings, conserving 
energy, momentum and angular momentum. Initially the A balls are placed 
randomly within a sphere of radius R — RqA^^"^, and the initial velocities 
are chosen as a superposition of thermal (Maxwell-Boltzmann) and collective 
motion. We use a spherically symmetric Hubble-like flow field for the initial 
collective motion: 

v(f} = -vof f/R (1) 

where wq/ is a model parameter, wq/ > for ingoing flow and vqj < for 
outgoing flow. We fix the total energy E = Efi + Eth, and vary the fraction 

rj of the flow energy, rj = Eji/E, where Eji = °i X]i=i ^ ^'^^ ^th are the 
flow energy and the thermal energy, respectively. Because of the way in which 
the system is built up, these energies will fluctuate from event to event with 
a relative uncertainty of the order of A~2. 

In our simulations we have chosen nuclear-scale parameters: m = 
940 MeV, The = 0.5 fm, Rq = 1.2 fm, < wq/ < 0.5 (in units of the ve- 
locity of light, c = 1), but since the behavior of the model only depends on 
the two combined parameters mv^t and rhc/ R, the choice of nuclear scale is 
not crucial. We choose A = 50, so with these parameters the initial radius of 
the system is 4.2 fm. 

We focus on four different types of event: • 'th20': The particles are 
started in 100% thermal motion inside a spherical container of radius 4.2 fm, 
at t = 20 fm/c the container walls are removed. • 'in': 100% ingoing 
flow. After interacting, the particles will move out again. This implosion- 
explosion process is intended to simulate some features of a heavy-ion col- 
lision. • '50/50': 50% thermal motion -f 50% outgoing flow, simulating an 
explosion from a non-thermalized state. • 'lOOout?': 100% outgoing flow in- 
side a spherical container of radius 7 fm. 

The results are averaged over an ensemble of 20 events of each kind. 

3 Thermodynamic considerations 

It is clear that we cannot use ordinary thermodynamics (or its well-known 
extensions to small systems □ or to small deviations from equilibrium cl) for 
the description of the overall behavior of our model. First, it is not clear 
that equilibrium prevails, even locally. Indeed we wish to investigate to what 
extent equilibrium is reached in the course of an implosion-explosion process. 
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Second, our system has no fixed volume, it expands freely into the vacuum. 
It is the combination of these two facts, no temperature and no volume- that 
makes our approach different from much previous work on the subject a. 

Equilibrium thermodynamics is linked to the motion of the individual 
constituents making up the macroscopic system via the entropy u'tl. A natural 
starting point for the investigation of the overall, i.e. the "thermodynamic", 
behavior of our system is therefore to apply an expression similar to the 
entropy, but in a way that makes sense in this highly non-equilibrium system. 

To study one-body observables, we reduce the QA dimensional phase-space 
of the A particles to 6 dimensions in the standard way □. Then we introduce 
a finite grid in the reduced phase-space, dividing each of the 6 axes into D 
segments. Instead of working with a fixed grid in phase-space, which would 
give us the usual entropy^, we let the entire grid expand or contract along 
with the swarm of points in phase-space in a uniform way: The outer grid 
edge follows the outermost point, the boxes are of equal size, and the number 
of boxes is kept fixed, thus the physical size {e.g. in units of ?i ) of each box 
in phase-space varies with time. This is to deal with the no volume problem, 
we mentioned above. We then introduce the pseudo-entropy as 

where 

number of points in box i 
total number of points in phase space 

and ^ is a normalization constant. We choose ^ as the theoretical maximum 
value of — X^i-Pj log(pi) so that E e [0, 1]. In the current set-up, for the case 
of N points in a reduced phase-space divided into D^ boxes, ^ is the smaller 
number of log(Af) and 61og(Z3). We refer to the quantity S as the pseudo- 
entropy instead of entropy, since some important features of the entropy, e.g. 
that it increases, are not retained in this formulation. Nevertheless, we shall 
see that S has some nice properties, including that of characterizing the degree 
of equilibrium. For a system of fixed volume in equilibrium, E is the usual 
one-body entropy (apart from normalization). 

4 Results 

In the calculations presented here, we have chosen I? = 7, so the 6 dimensional 
phase-space is divided into 7^ boxes. Because we are dealing with a small 



}_^P^l0g{p,) (2) 



"in the limit D ^ cxd in an ensemble of infinitely many events 
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number of events (typically 20) in the ensembles, the precise values of S 
depend on the choice of grid. We have, however, verified that S behaves 
qualitatively similar to what is shown here over a range of grid-sizes, varying 
D between 2 and 9. 

To give an idea of the dynamics and the timescales, we show in Fig. m the 
scattering rate. 




5 10 15 20 25 30 35 40 45 50 
time (fm/c) 

Figure 1. The scattering rate (number of scatterings per particle per time unit) for the four 
cases mentioned in the text. In the implosion-explosion event ('in') practically all particles 
scatter around the time t ~ 6 fm/c. This is the time when the system is maximally 
compressed. Then, as the expansion begins, the scattering rate decreases until t ~ 15 fm/c, 
when interactions have essentially ceased. In the 'lOOoutT' case the particles start to hit 
the container wall at i ~ 5 fm/c (the scatterings against the container walls are not counted 
here), and the peak in the scattering rate at i ~ 20 fm/c results from particles moving back 
after hitting the wall and scattering against other particles still on their way out. 

Fig. shows how the pseudo-entropy behaves in each of the four cases. 
In the 'th20' case, the pseudo-entropy S ~ 1 as long as the particles are in 
equilibrium at fixed volume inside the container. Then at t = 20 fm/c, when 
the container is removed and the system starts to expand, the pseudo-entropy 
decreases, reflecting the fact that the system goes out of equilibriuirn. 

In the case TOOoutT', where the particles are started in an extreme non- 
equilibrium situation, the pseudo-entropy is low (E = 0.8 is a low value in 
this context), but increases towards S = 1 as the scatterings equilibrize the 
system. By comparison with Fig. n^ one can see that the first jump in S at 



'' The particles stay almost thermalized, though, in the sense that they retain their Maxwell- 
Boltzmann velocity distribution. But they are certainly not in equilibrium, since this means 
that the phase-space distribution is independent of time. 
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Figure 2. The pseudo-entropy for the four cases described in the text. This variable seems 
to quantify the degree of equilibrium in the system, S = 1 characterising an equilibrium 
state. 

t ~ 5 fm/c is due to particles scattering against the container wall (when 
particles hit the wall their velocity is reversed, so in this process many new 
states in phase-space are being populated), and the second "jump" around 
t ~ 20 fm/c is due to the many particle-particle scatterings around this time. 

The interesting case is the implosion-explosion ('in') scenario, since here 
we do not know in advance if the system reaches a state of equilibrium or 
not. From the fact that the pseudo-entropy in Fig. reaches a value E ~ 1, 
the same value as the known equilibrium case 'th20', we infer that the system 
is in a state of equilibrium around i ~ 6 fm/c (which is also the time of 
m.aximum compression). We have checked that the speed distribution of the 
particles becomes nearly Maxwellian from i ~ 6 fm/c with a temperature 
in the compressed state of 47 MeV, which is also the theoretical value of 
the temperature assuming that all of the initial flow energy is converted to 
thermal energy. 

Another interesting feature of the pseudo-entropy is that it seems to decay 
in a characteristic fashion when the system expands from a state of equilib- 
rium, see Fig. y. 

5 Conclusions 

In this note we address the problem of thermodynamic equilibration in the 
context of heavy- ion collisions. We have defined a variable inspired by the 
entropy which, at least for the cases we have considered, seems to characterize 
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Figure 3. The pseudo-entropy in a log-log plot, together with the functions: 1.52 1~^'^ 
and 1.87 1""'^. The decrease of S during the initial expansion of the system seems to 
follow a power law when a state of equilibrium was present, in contrast to the case '50/50' 
(intended to simulate the expansion from a not-fuUy thermalized state) which does not 
show this behavior. At later times S decreases less rapidly and turns over to approach a 
finite limiting value, one sees the beginning of this behavior at the curve 'in'. 

the degree of equilibrium in an a priori highly non-equilibrium process such 
as an explosion. Now, more theoretical work needs to be done in order to 
understand why S behaves in this seemingly interesting way. 
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